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Strongly correlated electron liquids which occur in quantizing magnetic fields reveal 
a cornucopia of fascinating quantum phenomena such as fractionally charged quasi- 
particles, anyonic statistics, topological order, and many others. Probing these effects 
in GaAs-based systems, where electron interactions are relatively weak, requires sub- 
kelvin temperatures and record-high electron mobilities, rendering some of the most 
interesting states too fragile and difficult to access. This prompted a quest for new high- 
mobility systems with stronger electron interactions. Recently, fractional-quantized 
Hall effect was observed in suspended graphene (SG), a free-standing monolayer of 
carbon, where it was found to persist up to T = 10 K. The best results in those ex- 
periments were obtained on micron-size flakes, on which only two-terminal transport 
measurements could be performed. Here we pose and solve the problem of extracting 
transport coefficients of a fractional quantum Hall state from the two-terminal conduc- 
tance. We develop a method, based on the conformal invariance of two-dimensional 
magnetotransport, and illustrate its use by analyzing the measurements on SG. From 
the temperature dependence of longitudinal conductivity, extracted from the measured 
two-terminal conductance, we estimate the energy gap of quasiparticle excitations in 
the fractional-quantized v = 1/3 state. The gap is found to be significantly larger than 
in GaAs-based structures, signaling much stronger electron interactions in suspended 
graphene. Our approach provides a new tool for the studies of quantum transport in 
suspended graphene and other nanoscale systems. 



INTRODUCTION 

Fractional quantum Hall effect is a remarkable man- 
ifestation of electron interactions in two spatial dimen- 
sions 0, Q ■ The continued interest in this phenomenon 
is due to the large variety of fractionally quantized states 
and their rich physical properties Q. In particular, quasi- 
particles of such states can carry a fraction of electron 
charge 0, 0, [H] , and obey anyonic statistics, rather than 
the usual bosonic or fermionic statistics jy, 0]. It was 
predicted that non-abelian excitations may exist in the 
so-called v = 5/2 state @. So far most of experimen- 
tal studies of these effects have been conducted in high- 
mobility GaAs-based structures, where the electron inter- 
actions are relatively weak. This limits the temperatures 
at which the fractional Hall effect can be observed to 
T < 1 K. Many of the states of interest, in particular the 
5/2 state, were found to be extremely fragile, which is re- 
flected in the small energy gaps of elementary excitations. 
As a result, such states are only found in certain ultra- 
high- mobility GaAs structures @, [l(| , and probing their 
properties remains a challenge. This stimulated search 
for new two-dimensional systems with stronger electron 
interactions, which would host a larger variety of strongly 
correlated electron states. 

Recently, much interest was generated by the real- 
ization of a new two-dimensional electron system in 
gra phene, a one atom thick layer of crystalline carbon 
111 |l2j. One of the most exciting phenomena discov- 



ered in this material is the anomalous integer-quantized 
Hall effect (QHE) [H [3]. Due to the massless Dirac 
character of the carrier dispersion, the cyclotron energy 
in graphene can be orders of magnitude larger than in 
systems with massive charge carriers in similar magnetic 
fields, reaching a few thousand kclvin in 10 tesla. As a 
result, the QHE in graphene persists up to room tem- 
perature 15j. Another interesting feature is that, being 



a semimetal, graphene hosts a family of Landau levels 
with particle/hole symmetry, resulting in a particle/hole- 
symmetric arrangement of QHE plateaus in the trans- 
verse conductance, a xy — ve 2 /h {v = ±2, ±6, ±10...). A 
period-four regularity in these filling factors reflects the 
four-fold spin and valley degeneracy of Landau levels. 

Graphene is a truly two-dimensional material, in which 
electron interactions are expected to be stronger than 
in other systems. This should make interaction-induced 
QHE phenomena, such as the quantum Hall ferromag- 
netism (lifting of the valley and spin degeneracies of 
Landau levels) and the fractional quan tum Hall effect 
(FQHE) easily observable 0, [T3, H, lH] . Indeed, several 
integer QHE states outside the sequence v — ±2, ±6... 



have been reported [20|, l2l|, l22|. These states, how- 



ever, could only be observed in very strong magnetic 
fields. The fragile character of these states was linked 
to substrate-related disorder which, by inducing spatial 
variation of electron density, suppresses the effects of in- 
teraction. 

A breakthrough, however, was achieved very recently 
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FIG. 1: Theoretical fit (red) of the measured two-terminal 
conductance (black). The longitudinal conductivity a xx 
(blue), used in the fit, has peaks at the plateau-to-plateau 
transitions and minima on the plateaus (vertical dashed lines 
and arrows, respectively). The best fit was obtained using 
the semicircle relation between a xy and a xx , and treating the 
effective aspect ratio L/W, as well as the peaks' positions 
and widths, as fitting parameters as described in the text. 
Values at the minima, v = z/*, are used to estimate a xx for 
incompressible QHE states. 



in suspended graphene (SG) [23j, [24|. Electron mobil- 
ity in SG, enhanced due to the absence of substrate, was 
found to approach the ballistic limit in micron-size flakes. 
Low disorder makes SG an ideal system for studying the 
effects of interaction. In recent measurements on sus- 
pended bilayer graphene, interaction-induced QHE states 
were observed in fields below 1 tesla [25[ . Most recently, 



fractional Hall effect was reported in sin gle-l ayer SG, be- 
coming visible in fields as low as 2 tesla 2(| 27 1. 



Best results in these experiments could be achieved 
with micron-size flakes, since in larger flakes the electron 
density becomes nonuniform due to SG sagging. How- 
ever, in small flakes the standard four-probe measure- 
ment needed for separate measurement of the longitudi- 
nal and Hall conductivitiy is difficult to perform, because 
voltage probes in such small samples can short-circuit 
the Hall voltage. Observing the QHE using the standard 
Hall-bar geometry remains a challenge. In contrast, two- 
terminal measurements were found to consistently yield 
good results 24, 25, 2^, 27 1. However, while pointing to 
the existence of QHE states, the results of a two-terminal 
measurement are not amenable to straightforward inter- 
pretation. Thus, no quantitative characteristics of the 
new QHE states were obtained in Refs. 26|, 27 1. 



DECONVOLVING THE TWO-TERMINAL 
CONDUCTANCE 

In previous studies of QHE the longitudinal and Hall 
conductivities a xx and a xy , obtained from Hall-bar mea- 
surements, have been the quantities of choice, since their 
properties are most directly linked to the underlying 
physics. The width of the plateaus in a xy and the depth 
of the corresponding minima in a xx provide key infor- 
mation on localized states and on the electron interac- 
tion strength. Activation-like temperature dependence 
of 

(J xx CeLH be used to extract the excitation gap, while 
deviations from simple activation behavior can reveal the 
nature of electron transport mechanism. 

Since at present there are no reliable Hall-bar measure- 
ments in SG, it is tempting to use the two-terminal con- 
ductance for extracting the components of the conductiv- 
ity tensor a xx and a xy . However, the two-terminal con- 
ductance depends simultaneously on a xx , a xy and sample 
geometry, and thus 'deconvolving' it requires an addi- 
tional input. It was pointed out in Ref. 30] that such 
input can be provided by the conformal invariance of 
the magnetotransport problem. In this approach, a xx 
and a xy are interpreted as a real and imaginary part of 



a complex number a = a x 



>xy, 



and thereupon the 



transport equations become conformally invariant. Ap- 
plied to a rectangular two-lead geometry, theory yields 
a specific dependence of the two-terminal conductance 
on <j xx , a xy and the sample aspect ratio L/W. Inter- 
estingly, because of the conformal invariance, the same 
dependence describes the two-terminal conductance for 
an arbitrary sample shape, whereby the 'effective aspect 
ratio' encodes the dependence on sample geometry. 

Drawing on these observations about the role of confor- 
mal invariance, here we present a method for extracting 
transport coefficients a xx and a xy from the two-terminal 
measurements, and illustrate it using the data obtained 
as described in Ref. 26]. The measurements were car- 
ried out on suspended graphene samples for tempera- 
tures ranging from 1 K to 80 K and fields up to 12 Tesla. 
The samples were fabricated from conventional devices 
mechanically exfoliated onto Si/Si02 substrates by re- 
moving the Si02 layer with chemical etching [24| . In the 
final device the graphene sample is suspended from two 
Au/Ti pads which are split into two pairs one to apply 
the current and the other for probing voltage. The sam- 
ples were typically 0.6 /mi long and 1 — 3/im wide and 
were probed in a regime where transport was ballistic. 

The 'half-integer' QHE in these samples is observed at 
fields as low as IT. The electron interaction effects be- 
come important at fields above 4T, resulting in new quan- 
tized plateaus which correspond to integer QHE states 
at v = 0,±1,±3, and a FQHE state at v = 1/3 
The plateaus in conductance at these values of v ex- 
hibit nearly perfect quantization expected for the two- 
terminal conductance in the QHE regime at fields above 
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8T. At moderate fields (4 — 8T) the plateaus are less 
well developed, which is likely to be a result of spatial 
density inhomogeneity. The data obtained at the highest 
field, B = 1 2T, in which both the odd-integer and FQHE 
states are well developed, are best suited for our analysis. 

As illustrated in Fig[TJ our approach provides a good 
description of the data in a fairly wide range of filling 
factors. Here we focus on the v = 1/3 FQHE state and 
the v = 1 interaction-induced QHE state, and extract 
a xx on the plateaus. The temperature dependence of a xx , 
analyzed in terms of the activation transport mechanism, 
is used to estimate the excitation gap in these states. 
We find that the gap values which are two to five times 
greater than that measured in GaAs-based structures 1, 

Our approach, discussed in detail below, works best 
at temperatures which are neither too high, nor too low. 
The limitation from a high-temperature side stems from 
the nature of the semicircle relation, employed to describe 
the density dependence of a xy and a xx (see 33|, 34 1 and 
references therein) . This relation, which is also rooted in 
the complex- variable interpretation of magnetotransport, 
has been thoroughly tested in GaAs-based quantum Hall 



where the quantity I(b, a) is defined as an integral 



35, 36 



devices, both for integer and fractional QHE 
It was found to work well deep in the QHE regime, i.e. 
at not too high temperatures, and less well at elevated 
temperatures. 

At the lowest temperatures, the conductance of our 
SG samples exhibits pronounced mesoscopic fluctuations. 
These fluctuations dominate on the QHE plateaus, ren- 
dering the description in terms of an average <j xx inad- 
equate. Yet, as we show below, the approach based on 
the combination of conformal invariance and semicircle 
relation yields reasonably good results in a fairly wide 
range of intermediate temperatures, not too low and not 
too high. 

At temperatures such that the mesoscopic fluctuations 
are not too prominent, we adopt an approximation in 
which the transport coefficients a xx and a xy take con- 
stant values throughout the sample. For a sample of 
a rectangular shape with ideal contacts at the opposite 
sides of the rectangle the solution to the transport prob- 



lem has been long known |28|, [29[. The result can be 
summarized in a compact form following Ref . [28[ , where 
it was derived with the help of conformal mapping. For 
a rectangle of length L and width W it is convenient to 
parameterize the aspect ratio as 



where K denotes the complete elliptic integral of the first 
kind. The resistance of such a rectangle is then given by 



n(p xx , p xy ) — y p xx -t- p xy jt^ _^ , K — — 



>1. (2) 



7(6, a) 



i(e-i)(*'+oi a -i(£ -Oh* 

a 

a ± = l/2±e/n, (4) 

where 9 = a,Tcta,n(p xy / p xx ) is the Hall angle. Because 
of conformal invariance of the two-dimensional mag- 
netotransport problem, Eq.@ also describes the two- 
terminal resistance for a sample of an arbitrary shape, 
with the dependence of R on the sample geometry char- 
acterized by a single parameter £, the aspect ratio of an 
'equivalent rectangle' [3(ij . 

It is instructive to consider the behavior at large Hall 
angles, p xx <C \p xy \, 9 — ► ±7r/2, which is a regime rele- 
vant for our discussion below. Without loss of generality 
we can assume p xy > 0, in which case in the limit p xx — > 
we have a_ — > and a+ — » 1. In this limit, the integrals 
in Eq. ([2]), dominated by the regions near £ = — 1 in 
7(1, —1) and £ = k' in I(k', 1), diverge as l/a_. Extract- 
ing the leading contributions, we find that they are the 
same in both cases: 



'(1,-1) 



1 



(k' + 1)( 



/(fc',1) 



(k' + l)a- 



(5) 



Thus the ratio I(k', 1)/7(1, — 1) tends to one, giving 
72 = \Pxy\i which is the behavior expected in the dis- 
sipationless QHE state, when p xx = 0. By carrying out 
the expansion in small a_ to next order, a contribution 
linear in p xx can be found (see Ref. [28|): 



7? 



Pxx9{t), g{t) = h\ 



1 



2Vk 



(6) 



The quantity g(£) is positive for £ = L/W > 1, and neg- 
ative for £ = L/W < 1; it vanishes at I = 1. For I » 1 
(long and narrow sample) the function g{€) is approxi- 
mately linear: g(£) w £. For I 1 (short and wide sam- 
ple) the function g(£) behaves as l/£. This behavior is 
consistent with the results expected for uniform current 
flow, R(£ > 1) - p xx £, R(£ < 1) ~ 1/((t xx £). Deviation 
from a quantized conductance value, described by Eq.©, 
offers a way to extract a xx from two-terminal measure- 
ments. 

Given the values of a xx , <r xy and £, the resistance 7? can 
be obtained by numerical evaluation of the quantities in 
Eq.([2]). However, while the integral ([3|) converges for all 
—tt/2 < 9 < tt/2, the convergence is slow for large Hall 
angles, a xx -C \cxy\, 9 —> ±tt/2, because of the power- 
law singularities of the integrand. Because of that, we 
found it more convenient to evaluate 7? using another 
method which was developed by Rendell and Girvin [i^] . 
In this approach, the current density is found as an expo- 
nential of certain infinite sums. The total current 7 and 
source-drain voltage Vsd are obtained by integrating the 
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current density and electric field over appropriate con- 
tours, after which the resistance R is found as the ratio, 
R = Vsd/I- The results obtained by this method are 
identical to those found from Eq.([2]). However, since the 
infinite sums giving current density converge rapidly, the 
numerics turns out to be substantially simpler than when 
Eq.([2]) is used directly. In what follows, we will use the 
approach of Ref. [29( to evaluate the resistance. 

There are several ways to use this approach for de- 
termining <j xx . One is to focus on the plateaus, where 
c xx is small and a xy is quantized, a xy = ve 2 /h. Ex- 
panding G in the small ratio cr xx /\a xy \ <C 1, the devi- 
ation from a quantized value can be expressed as G = 
\v\e 2 /h — g(£)a xx + 0(cr xx ), where the coefficient g(£) is 
a function of the aspect ratio only [see Eq.©]. Despite 
a conceptual simplicity of this approach, we found it dif- 
ficult to implement, since the effective L/W value may 
significantly deviate from the geometric aspect ratio of 
the sample, and thus should be treated as a fitting pa- 
rameter [3l[ . Further, since a xx and a xy change with v, 
the conductance plateaus exhibit N-shaped distortions, 
rendering the deviation in G from a quantized value un- 
suitable for accurate estimation of a xx - 

Considerably more reliable results can be obtained by 
focusing on the N-shaped distortions of QHE plateaus, 
since matching the theoretical model to an entire curve 
G{v) puts a substantially more stringent constraint on 
the fitting parameters. The N-shaped features can be 
described by the density-dependent a X v(v) and a X x{v), 
which obey the semicircle relation 3(1, 31| . For a plateau- 
to-plateau transition between incompressible filling fac- 
tors V\ < V2 this relation gives a 2 x = {a xy — vi)(v2 — a xy ) 
(in units of e 2 /h). For fitting the conductance data 
shown in Fig[TJ which exhibits incompressible states at 
v — 0, 1/3, 1, 2, we model the contribution of each QHE 
transition by a Gaussian, 



V X <V C < l/ 2 , 



(7) 



and find the corresponding a xy from the semicircle rela- 
tion. This gives a contribution to the longitudinal and 
Hall conductivity of each of the relevant Landau levels 
or sublevels. The net conductivity a xy (y) and cr xx (y), 
found as a sum of such independent contributions (blue 
curve in FigO}, is then used to calculate the dependence 
G{v) = 1/R, where R is given by Eq.©. 

We treat the a xx peak positions and widths, as well 
as the effective ratio L/W, as variational parameters. As 
illustrated in Fig[TJ the Gaussian model with individually 
varying peak widths and positions provides a rather good 
description of the data. The sum of Gaussian peaks gives 
the quantity a xx (i>). The approach based on treating 
L/W as a variational parameter, in general different from 
the actual sample aspect ratio, works rather well in the 
integer QHE regime 3l|. It was conjectured 3l| that 
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FIG. 2: The statistical procedure used for extracting a xx is 
illustrated for v = 1/3 plateau of T = 6K conductance trace 
(see Fig|3h). Scanning a range of fitting parameters, includ- 
ing the aspect ratio L/W and the widths and positions of 
gaussian peaks, Eq.([7|, yields several thousands fits. Quality 
of the fits is characterized by the standard deviation, Eq.©. 
The value at the minimum, v = f„, (see Fig.([T} gives con- 
ductivity (Txx for each fit. The resulting conductivity was 
estimated by averaging the values of a xx over the group of 
fits within 30% of the best fit, as shown by a box in the in- 
set. Error bar was determined from the standard deviation 
of a xx within this group of fits. Similar results can be found 
from the quantity Ave„(a xx ), obtained by averaging a xx over 
n best fits, extrapolated to n = 1. 



sample-dependent specifics of the current flow pattern 
such as those due to imperfect contacts and/or contact 
doping. 



STATISTICAL ANALYSIS OF THE FITTING 
PROCEDURE 

In order to extract a xx in the incompressible v = 1/3 
and v — 1 states, we analyzed a set of fits which best 
follow the data near the corresponding plateaus (FigQJ. 
In both cases we found an optimal effective aspect ratio 
L/W w 0.59 (such deviation from the geometric aspect 
ratio, which in the sample [26| was close to 0.4, is con- 
sistent with the results of Ref.[31(). Best fits were found 
from the square-deviation in the conductance averaged 
over a range of densities on the plateau and around it, 



(8) 



variations in the best-fit value of L/W account for the 



In the case of the 1 /3 plateau, a small interval near v w 
0.4, which is probably related to another incipient FQHE 



feature, was excluded from the integral in Eq.©. The 
values at the minima, v = (see FigJTJ) , were taken 
as an estimate of the longitudinal conductivity a xx of 
incompressible QHE states. 

Search for the best fit was performed by optimiz- 
ing fitting parameters, which include the positions and 
widths of Gaussian peaks, Eq.®, and the aspect ra- 
tio L/W. Our statistical analysis, based on comparing 
square-deviations D, Eq.©, for several thousands differ- 
ent, randomly chosen parameter values, is illustrated in 
FigfS] The value of a xx was estimated by averaging over 
a group of fits within 30% of the best fit, marked by a 
box in Figf2] inset. Statistical error was estimated from 
the spread in the a xx values found from the fits within 
this group. The resulting error bar, displayed in Fig|3]b, 
is below 10% at 10K, and increases to 20-25% at 1.2K 
as a result of pronounced mesoscopic fluctuations devel- 
oping at low temperatures. Alternatively, the value of 
a xx could be determined be analyzing the mean value 
Ave n (a xx ) taken over n best fits, ordered according to 
their square-deviation D, as shown in FigO Linear ex- 
trapolation of Ave n (a xx ) to n = 1 yields results which 
are close to those found by averaging over the group of 
best fits, as described above, and has an advantage of 
being less subject to statistical error. 

RESULTS AND DISCUSSION 

To extract a xx at different temperatures, we used the 
12 Tesla, displayed in FigEK- On the plateaus, the con- 
ductance increases with temperature faster than off the 
plateaus, consistent with the behavior expected in the 
QHE regime. Both plateau features remain visible up 
to T ~ 20 K. For the highest temperatures, however, 
the density-dependent conductivity <J xx {v) exhibits very 
shallow minima on the plateaus, giving rise to a large 
statistical error in fitting. We attribute this behavior to 
the temperatures T > 20 K lying outside the range of 
applicability of the semicircle model, similar to the high- 
est temperatures in Refs.Jil, |36j]. Thus we exclude the 
T = 20, 40 K traces from the analysis. Statistical er- 
ror, estimated from fluctuations of a xx over a group of 
best fits, as discussed above, gives the error bars shown 
in Figj3)3. Larger statistical error found for T = 1.2 K 
reflects growth of mesoscopic fluctuations, which domi- 
nate the transport on the QHE plateaus at the lowest 
temperatures. 

The conductivity values on the plateaus, obtained from 
the conductance traces with T = 1.2, 4.2, 6 and 10 
K, were analyzed in terms of the activation behavior, 
a xx ~ exp(-A/2fc B T) (see FigSp). The best fit val- 
ues of the energy gap are A/ks = 10.4 K for v = 1 
and A/k B = 4.4 K for v = 1/3. The relatively high 
energy scale is consistent with the disappearance of the 
corresponding plateaus at about 20K. We also considered 
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FIG. 3: a) Temperature dependence of the two-terminal 
conductance from which a xx on the v — 1/3 and v = 1 
plateaus is extracted, b) Fits of the temperature dependence 
of Oxx to the activation model, a xx ~ exp(— A/2ksT). The 
best-fit values of the activation gap are A„=i = 10.4 K and 
A„ = i/3 = 4.4 K. The procedure used to estimate error bars 
is outlined in Fig[2] 



fits to the variable-range-hopping dependence, a xx ~ 
exp (— (T^/T) 11 ), j] — 1/2, but have not found a dis- 
cernible statistical advantage over the activation depen- 
dence. 

It is instructive to compare these results with the ac- 
tivation gap measured in the v — 1/3 state in GaAs 
structures. In the first experiment, in which the FQHE 
at v = 1/3 was discovered 1], the temperature depen- 
dence of a xx was weak, and the activation gap could not 
be determined. Subsequent measurements on samples of 
higher mobility [12] revealed activated behavior of con- 
ductivity at temperatures T > , with of the order of 
0.1 — IK, depending on the field. The activated depen- 
dence crossed over to the variable-range-hopping behav- 
ior at lower temperatures, indicative of the localization 
of quasiparticles. Despite very high mobility of GaAs 
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structures, the extracted gap value, which was about 2 K 
at B = 12 T, is about 2.5 times lower than the value for 
SG obtained above, indicating a more robust nature of 
FQHE in graphene. 

Given that the FQHE features in SG remain clearly 
visible up to 10-20 K (see Fig[3h.), whereas in the mea- 
surements on GaAs the FQHE persisted only up to about 
2 K, one may expect subsequent measurements on cleaner 
SG samples to revise the gap values. Another indi- 
cation that our estimate of the gap is merely a lower 
bound arises from comparison to theoretically predicted 
values. Theoretical estimate for v = 1/3 in SG gives 
A w=1 /3 = ae 2 /n£B, where a s» 0.1 [17]. taking the di- 
electric constant K — 5.24, which is the RPA result for 



intrinsic screening function of graphene [37|, for B = 12 T 



we obtain A,, 



= 1/3 



42 K. 



This situation can be contrasted with GaAs, where the 
theoretical value, A„ =1 / 3 w ae 2 /eIb, with a s» 0.03 and 
e = 12.8, is only a factor of 2 .5 g reater than experimental 
value measured at B = 12 T 32J. A relatively small value 
of the prefactor a accounts for the effect of finite width of 
GaAs quantum wells [IH, which makes the short-range 
interactions weaker than in a truly two-dimensional sys- 
tem such as graphene. In higher fields, B ~ 20 T, the 
theoretical limit in GaAs has been nearly reached [32j ]. 
due to the enhanced role of interactions compared to dis- 
order. 

We therefore believe that the departure of the gap in- 
ferred in present work from theoretical predictions re- 
flects the effect of disorder present in the system. In 
particular, it was pointed out that rippling of a sus- 
pended graphene sheet may result in formation of local- 
ized midgap states 12j. Yet, similar to the case of GaAs, 
the effect of disorder should become weaker at higher 
magnetic fields. This would bring the FQHE gaps closer 
to the very large theoretically predicted values, roughly a 
factor of 20 greater than those in GaAs. Realizing FQHE 
with larger gaps should also be possible at lower magnetic 
fields once the sample quality is improved. 



SUMMARY 

In this work, we developed a general method for ex- 
tracting transport coefficients in a QHE state from the 
two-terminal conductance. We demonstrate that, while 
the two-terminal conductance depends in a fairly com- 
plicated way on the sample geometry, as well as on a xx 
and <r X u, a reliable procedure for determining these values 
can be developed. This type of analysis is made possible 
by constraints on this quantity arising from conformal 
invariance of the two-dimensional magnetotransport and 
the semicircle relation between a xx and a xy in a QHE 
state. We apply our approach to analyze the fractional 
and integer QHE states in suspended graphene flakes, 
where, because of small sample size, only two-terminal 



measurements can be performed 26j, |27j ■ 

We estimate the energy gap of the quasiparticle ex- 
citations in the v = 1/3 FQHE state by analyzing the 
temperature dependence of a xx . The gap is found to be 
significantly larger than in semiconducting systems, sig- 
naling stronger electron interactions in graphene. From a 
comparison to measurements in GaAs and to theoretical 
estimates, we conclude that the effects of electron inter- 
actions in current SG samples, despite being somewhat 
masked by disorder, are stronger than in high-quality 
GaAs structures. We expect that the future experiments 
on cleaner SG will reveal exceptionally robust FQHE, 
with gaps reaching a few tens of Kelvin. Given the rich- 
ness and diversity of FQHE phenomena in GaAs struc- 
tures, it is natural to expect that SG will be shown to 
host new types of FQHE, which are not observable in 
semiconducting structures because of the much weaker 
Coulomb interactions. 
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